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Abstract 

Aims. This paper is the first of a series investigating Lyman-alpha (hereafter Lya) radiation transfer through hydrodynamical simula- 
tions of galaxy formation. Its aim is to assess the impact of the interstellar medium (ISM) physics on Lya radiation transfer and to 
quantify how galaxy orientation with respect to the line of sight alters observational signatures. 

Methods. We compare the results of Lya radiation transfer calculations through the ISM of a couple of idealized galaxy simulations 
in a dark matter halo of ~ 10 10 M o . In the first one, Gl, this ISM is modeled using physics typical of large scale cosmological hydro- 
dynamics simulations of galaxy formation, where gas is prevented from radiatively cooling below 10 4 K. In the second one, G2, gas 
is allowed to radiate away more of its internal energy via metal lines and consequently fragments into dense star-forming clouds. 
Results. First, as expected, the small-scale structuration of the ISM plays a determinant role in shaping a galaxy's Lya properties. 
The artificially warm, and hence smooth, ISM of Gl yields an escape fraction of ~ 50% at the Lya line center, and produces sym- 
metrical double-peak profiles. On the contrary, in G2, most young stars are embedded in thick star-forming clouds, and the result 
is a ~ 10 times lower escape fraction. G2 also displays a stronger outflowing velocity field, which favors the escape of red- shifted 
photons, resulting in an asymmetric Lya line. Second, the Lya properties of G2 strongly depend on the inclination at which it is 
observed: From edge-on to face-on, the line goes from a double-peak profile with an equivalent width of ~ -5 A to a 15 times more 
luminous red- shifted asymmetric line with EW ~ 90A. 

Conclusions. The remarkable discrepancy in the Lya properties we derived for two ISM models raises a fundamental issue. In effect, 
it demonstrates that Lya radiation transfer calculations can only lead to realistic properties in simulations where galaxies are resolved 
into giant molecular clouds. Such a stringent requirement translates into severe constraints both in terms of ISM physics modeling 
and numerical resolution, putting these calculations out of reach of current large scale cosmological simulations. Finally, we find 
inclination effects to be much stronger for Lya photons than for continuum radiation. This could potentially introduce severe biases 
in the selection function of narrow-band Lya emitter surveys, and in their interpretation, and we predict these surveys could indeed 
miss a significant fraction of the high-z galaxy population. 
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1. Introduction 

In the last decade, the Lyman-alpha (Lya) emission line has 
become an observational tool of choice to detect high redshift 



galaxies via narrow-band surveys ( e.g.|Hu et al.|19 98 ; Kudritzki 
et al.pOOQUShimasaku et al.|2006l|Ouchi et al.|2008|l2010[|Hu 
et al. |2010|) or blind spectroscopic searches (e.g. [v an Breukelen 
et al.|2005| |Rauch et al.[2008l|Cassata et al.|[2011| ). Today, the 
number of galaxies detected in this fashion (hereafter Lyman- 
Alpha Emitters, LAE) is becoming statistically significant, and 
LAEs play a major role in our census of high-z galaxies. At 
the same time, spectroscopic follow-ups of UV-selected galax- 
ies shed more and more light on the physical nature of LAE and 
on their place in the cosmic history of galaxy formation ( Shapley 
et al.|2003[ |Tapken et al.|2007| |Bielby et al.|20TT] ». One of the 
major challenges in the years to come, both theoretical and ob- 



servational, is yet to understand the details of the Lya line pro- 
files we observe: How do they relate (if they do) to any physical 
property of high-z galaxies ? 

Although a number of semi-analytic models for Lyman- 
alpha Emitting galaxies (LAEs) have been published (e.g. [Le| 
Dell iou et al.||2Q06| [QrsTet al.||2008[ [Dayal et al.||2008| [2009f 



Qrsi et al.|2011HDayal et al.|2011[|Garel et al.|2012| ), the com- 
plete radiation transfer through the interstellar medium of Lya 
emitting galaxies h as been taken into account in only a handful 
of previous studies fl¥s itsiomi 2006 ; Laursen et al. 2009 ; Barnes 



et al.|20ir||Yajima et al.|2011a|b see Table [1] for a summary.). 



In Tasitsiomi (2006), Lya radiation transfer is post- 



processed in the brightest Lya emitte r of a gas -dynamic s+N- 
body adaptive refinement tree (ART, [Kravtsov an d Gnedin 2005 ) 
simulation at z ~ 8, in order to investigate its detectability. In 
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Laursen et al.| ( |2009] ), a sample of nine galaxies at z = 3.6 taken 
from cosmological N-body/hydrodynamical TreeSPH simula- 
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cosmo zoom 


cosmo zoom 


cosmo zoom 


cosmo zoom 


nb of objects 


2 


1 


9 


3 


~ 950 


stellar mass 


1.8 x 10 9 M o for Gl 


~ 10 10 M o 


6 x 10 6 to 


1.5 x 10 10 M o 


4.3 x 10 9 M o 




4.9 x 10 8 M o for G2 




3 x 10 10 M o 


1.5 x 10 u M o 


9.3 x 10 9 M o 










7.5 x 10 n M o 


4.1 x 10 10 M o 


stellar mass resolution 


1.4 x 10 3 M o for G2 


2 x 10 4 M o 


1O 6 M 


not available 


1.9 x 1O 4 M 




7.7 x 10 3 M o for Gl 










spatial resolution (a) 


1 8 pc for G2 


29 pc 


137 pc w 


514 pc 


342 pc 




147 pc for Gl 










Hi temperature 


10 2 to 10 5 K 


10 3 to 10 4 K 


10 4 K 


10 43 to 10 5 K 





Table 1. Comparison of the 4 published studies of Lya radiation transfer through the interstellar medium of galaxies with this study. 

Resolution, in physical pc. This is either the minimum cell size, for AMR codes, or the gas gravitational softening length for 
SPH codes. In both cases, this reflects the smallest scale onto which a gas overdensity may feel its own gravity. 
^ This resolution corresponds to their S87 simulation, which best matches our halo mass. 



tions ( Sommer-Larsen 2006) are considered, sampled in mass. 
The main result of the paper is an anti-correlation between Lya 
escap e fraction and the mass of the galaxy. In Bar nes et al. 
(2011), three halos at z = 3 are selected in cosmological hy- 
drodynamic simulations (GADGET-2) aimed at reproducin g the 
physical p ropertie s of the host galaxie s of DLAs atz ~ 3 ( |Tescari| 
et al.|2009| ). In |Yajima etaL]p01 la|b] ), the same halo is followed 
at different redshifts, in order to study the evolution of the Lya 
properties of their galaxies with time. 

These previous studies are done through a warm interstellar 
medium in which the gas is cooled down to T~ 10 4 K, with the 
consequence that the formation of small scale structures is not 
modelled. Indeed, the pressure support of this warm gas prevents 
it from collapsing at scales smaller than its Jeans length, explain- 
ing the difference in spatial resolution between the different ex- 
periments (see section [2] for more details). However, theoreti- 
cal expectations suggest a strong dependance of the Lya trans- 
fer on the structure, and geometry of the interstellar medium of 
galaxies, and the main goal of our study is to investigate this 
point. Furthermore, two studies are dust-free, which prevent 
them from studying the Lya escape fraction from their config- 
urations. Finally, monochromatic approaches of the problem 
do not allow to derive Lya equivalent widths, and compare the 
transfer of continuum versus line photons. 

To overcome these limitations, we post-process hydrody- 



namical simulations of galaxy formation described in Dubois 
and Teyssier] ( |2008| ) perfor med with the RAMSE S code 



( |Teyssier||2002| ), with McLya ( Verham me et al.||2006| ), includ- 
ing Lya+continuum radiation transfer in a dusty medium, for 
two different ISM models : 1/ the reference model Gl, compa- 
rable to previous studies, where the gas is cooled down to 10 4 K, 
2/ a more realistic ISM model G2, where the gas is allowed to 
cool down to 100K, and the formation of small scale structures 
is followed. 

The plan of this paper is the following. We start to describe 
the hydrodynamical simulations used to post-process Lya radi- 
ation transfer. Then, we describe the radiative transfer of Lya 
photons in the hydrodynamical simulations with McLya. The 
fourth part presents a comparison of the Lya properties for the 
two ISM models. In the fifth part, we discuss the effect of ori- 



entation on the Lya properties of G2, which appears as an in- 
teresting bonus result of this work. The sixth part describes the 
diffuse Lya halo around G2. The last part summarizes the main 
conclusions. 



2. Description of the hydrodynamical simulations 

The results presented in this paper are based on the analysis of 
a couple of idealized high-resolution hydrodynamical simula- 
tions which follow the formation and evolution of an isolated 
disc galaxy embedded in a live dark matter (DM) halo. Both 
simulations were run with RAMSES (Teyssier 2002), using sub- 
grid physics modules described in Rasera and Teyssier (2006) 
and |Dubois and Teyssier| ( [2008] ). 



2.1. Initial Conditions 

We choose to follow the formation of a galaxy embedded in a 
halo of total mass (dark matter plus gas) M200 = 10 10 M o . The 
choice of this rather small halo mass is motivated by the expec- 
tation that most small-mass haloes host LAEs, while only a frac- 
tion of m assive galaxies emit Lyman-alpha at all (see e.g. |Garel| 
|et al.|2012l and references therein). 

We assume a NFW density profile ( Nava rro et al.|1996| ) with 
concentration parameter c = 10 for the DM halo. The 'virial' 
radius 7? 2 oo ~ 35 kpc of this latter is defined so that it encloses 
an average density equal to 200 times the mean matter density 
of the Universe at z = 0, assuming Q m = 1 and h = 1 . N ote that 
with the currently favored cosmology (e.g. Komatsu et al.|201l] ), 
this value of 7? 2 oo corresponds to a z ~ 1 halo of the same mass. 

We follow [Dubois an d Teyssier ([2008) to generate our ide- 
alized initial conditions and sample the DM halo with particles 
instead of using a static gravitational potential. This is necessary 
to allow gas and stars to exchange angular momentum with the 
DM. This process is particularly important in simulations where 
gas fragmentation occurs in the disc because the resulting star 
forming and gravitationally bound clouds should be driven to 
the bottom of the potential well by dynamical friction. In prac- 
tice, we sample our halo out to a radius of 3.2 x 7? 2 oo with 10 5 
DM particles of mass -1.8 1O 5 M . 
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Gl 


G2 


Stellar mass M* 


2.1 x 10 y M o 


6.8 x 10*M o 


Gas mass M gas 


5.0 x 10 8 M o 


3.0 x 10 8 M o 


M*/M gas 


~4 


~ 2 


Dust mass M dust 


8.7 x 10 6 M o 


6.3 x 10 6 M o 


SFR 


O.3M /yr 


l.lM /yr 



Table 2. Comparison of the physical properties of Gl and 
G2. From the SFR, and assuming a Salpeter IMF, we can de- 
rive intrinsic Lya luminosities of our galaxies: L(Lya)Qi = 
3.3 x lO 41 ^.^ 1 andL(Lya) G2 = 1.2 x 10 42 erg. s~ l . 



In order to form a centrifugally supported gas disc, we give 
this distribution of DM particles a specific angular momentum 
profile j(r) = y max x M(< r)/M 2 oo> where r is the distance to the 
halo center. The normalization y max is fixed by setting th e dimen- 
sionless s pin parameter of the halo to A = 0.04 (see e.g. Bullock 
et al.|200l) . The initial gas density profile is obtained by scaling 
that of the DM by the universal baryonic fraction fi> = 0.15. The 
gas velocity field is chosen to be identical to that of the DM. We 
enforce hydrostatic equilibrium, which uniquely defines the gas 
initial temperature (or pressure) profile, and its chemical com- 
position is pristine (76 % Hydrogen and 24 % Helium in mass). 



The simulated volume is a cubic box of size Lb ox = 300 
kpc on a side, sufficiently large to self-consistently follow the 
collapse of the gaseous halo for up to 6 Gyrs. Open bound- 
ary conditions let the galactic wind stream out of the box freely. 
The snapshots are chosen at the time where the galactic wind is 
reaching the virial radius, i.e. after 3 Gyr for Gl and 1 Gyr for 
G2. Although ages are different, M gas and Md us t are compara- 
ble for Gl and G2 (See Table [2]). The adaptive mesh refinement 
(AMR) grid is composed of 64" cells on the coarsest level, which 
corresponds to a cell size of ~ 4.7 kpc. Levels of refinement are 
triggered when at least one of the two following criteria is ful- 
filled: 

1. The total baryonic mass in a cell becomes larger than 
6 10 4 M Q , or a cell contains at least 8 DM particles 

2. A cell size becomes larger than 0.25 A] ( [Truelove et al.|1997] ) 

Refinement is allowed up to a maximum level which depends on 
the s imul ation. In our Gl (resp. G2) simulation (see Tab. [T]and 
Sec. 2.2), the maximum level allowed is 11 (resp. 14), which 
corresponds to a spatial resolution of 147 pc (resp. 18 pc). 



2.2. Physics of galaxy formation 

In this subsection, we briefly discuss how the physics (be- 
sides gravity and hydrodynamics) relevant to galaxy formation 
is modeled in the simulations. Table [2] summarizes the global 
properties of Gl and G2. 

Cooling: Radiative energy losses, assuming gas is in collisional 
ionization equilibrium, are computed in each grid cell using the 
metallicity depend ent cooling functions tabulated in Suthe rland] 
and Dopita ( |l993| . This allows the gas to cool down to 10 4 K in 
runs Gl and G2. In G2, we also account for the extra loss of en- 
er gy provided by metal line c ooling, following the prescription 
of Rosen and Bregman]( 1995 ). This allows the temperature of 
the gas to drop further down, to 100K in this run. 



of state (EoS) in high-density regions 

\Po 



(1) 



where p is the polytropic index, po is the gas density above 
which the polytropic law applies, and To is the minimum tem- 
perature of the gas due to cooling. We adopt a p = 2 polytropic 
index that ensures that the Jeans length is independent of the gas 
density for densities above po. For run Gl, To = 10 4 K and 
Po = 0.1 H.cm -3 , and the Jeans length of the gas is Aj ^ 5.8 kpc, 
comparable to the size of the entire galactic disc. This means 
that the ISM gas does not fragment, but instead is smoothly dis- 
tributed throughout the disc (see Fig[T]). For run G2, T Q = 100 K, 
and po = 10 H.cm -3 , so the constant Jeans length given by the 
polytropic EoS is Aj = 58 pc, which we resolve with ~ 3 cells of 
size Ax = 18 pc. This choice of T and po creates a multiphase 
medium with cold and dense bounded regions of size ~ 100 pc 
similar to the large giant molecular clouds (GMCs) of the Milky 
Way, and a low-density warm medium with strong turbulent mo- 
tions (see Fig [2]). 

Star for mation: Star formation is modeled as in IRasera andl 
Teyssier| ( |2006| ), with a random Poisson process spawning star 



cluster particles according to a classic Schmidt law. In other 
words, the star formation rate p* scales with the local gas density 
p as 

P* = e*f, (2) 

where % is the free fall time and 6* is the star formation effi- 
ciency. We adopt 6* = 1% in runs Gl and G2, which is in fair 
agreement with obse rvations over the gas density range we span 
(Krumholz and Tan 2007 ). Star formation is only permitted to 
occur in grid cells where the gas density is larger than p . Each 
star-cluster particle created is given a mass which is an integer 
multiple of the minimal mass m* = po x Ax 3 , where Ax is the 
size of a cell on the highest level of refinement. Therefore, our 
stellar mass resolution is m* ^ 1.4 10 3 M© for the G2 simulation 
(and m* 7.7 10 3 M© for the Gl simulation). We forbid the star 
formation process to consume more than 90% of the gas mass in 
the cell where it takes place during a time step. 



Feedbac k: We mo del supernovae (SN) feedback as in |Dubois 
and TeyssTer|(|2008]), i.e. we deposit kinetic energy in a sphere 



centered on the explosion, along with mass and momentum dis- 
tributed according to a Sedov-Taylor blast-wave profile. This 
kinetic approach ensures the formation of large-scale galactic 
winds in low-mass haloes, and galactic fountains in the most 
massive ones even at low resolution ( Dubois and Teyssier|2008 



Interstellar Medium: In order to prevent numerical fragmenta- 
tion, we assume the ideal gas transitions to a polytropic equation 



Such a precaution is necessary since it is w ell known that spu- 
rious thermal energy losses are catastrophic ([Navarro an d White 
|1993| ). We assume a standard Salpeter (1955) Initial Mass 
Function (IMF), with t] S n = 10% (in mass) of the newly formed 
stars finishing their ~ 10 Myr lives as type II SN. We further 
assume that 10 51 ergs of energy are released by each individual 
SN explosion with a typical ejected mass of 10 M©. SN are also 
responsible for releasing metals into their surroundings with a 
constant yield y = 0. 1 . These metals are passively advected. 



3. Lya transfer 

We use an improved version of the Monte Carlo radiation trans- 
fer code McLya of Verhamm e et al.| ([2006) to post-process the 
two simulations Gl and G2 described in Sec. |2] The most no- 
table technical improvement over the original McLya is that it 
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Figure 1. Left : Gl run galaxy seen face-on. Right : Gl run galaxy seen edge-on. The upper panels show the distribution of the 
gas density (in log cm~ 3 units), and the bottom panels the distribution of stellar density for stars younger than 10 Myr (in arbitrary 
log units) . 



now fully exploits the AMR grid structure of RAMSES, which 
allows us to perform the radiative transfer at the same spatial res- 
olution than our hydrodynamics simulations. The new version 
of McLya also features more detailed physi cs of the Lya line and 
UV continuum transfer, which include (see Scha erer et al.|201 1) : 
angular redistribution functions taking quantum mechanical ef- 



fects for Lya into account (Dijkst ra and Loeb||2008[ |Stenflo 



1980| ), frequency changes of Lya photons due to the recoil ef- 



fect (e.g. Zheng and Miralda-Escude 2002 ), the presence of deu- 
terium (a ssuming a canonical abundance of D/H = 3 x 10" 5 , 
Dijk stra et al.||2006| ), and anisotropic dust scattering using the 
Henyey-Greenstein phase function (using parameters adopted in 
|Witt and Gordon|2000| ). 

Before McLya can be used to process our simulations, we 
have to extract a set of gas properties for each cell which are 
not directly predicted by RAMSES: the gas velocity dispersion 



due to temperature and small-scale turbulence (the macroscopic 
velocity and temperature of each cell are known), the ionization 
state (or equivalently the neutral Hydrogen density), and the dust 
content. We discuss how we compute these quantities in the 
following sub-sections, and then explain our strategy to sample 
Lya emission. 



3.1. lonisation state of the gas 

The simulation outputs provide us with the total density p of 
gas in each cell, including H, He and metals. We derive the 
numerical density of H as = X H p/mn, where X H - 0.76 
is the mass fraction of H, and m H the mass of an H atom. Lya 
photons however only scatter onto neutral Hydrogen atoms (Hi), 
the numerical density of which we compute as nui = xnn, where 
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Figure 2. Same as Fig. 0for the G2 run. 



the neutral fraction x is evaluated assuming collisional ionization 
equilibrium (CIE) and hence only depends on temperature T. 

Using the polytropic equation of state (Eq. [TJ to prevent 
artificial fragmentation in the dense ISM unfortunately makes 
the gas temperature in these regions unknown. In all polytropic 
cells, we therefore set the temperature of the gas to Tq (= 10 4 K 
in Gl and 10 2 K in G2). Since we are assuming CIE, these tem- 
peratures imply that all the polytropic gas is neutral. 



lence, one has to quadratically add the turbulent velocity (v tur b) 
contribution and write b = (v 2 h + v 2 urb ) 1/2 . 

In practice, we set v tur b = 10 km s" 1 in both Gl and G2 
simulations. This value is the mean value of turbulent velocities 



computed in simulations of SN-driven turbulence by |Dib et al. 
( 2006 ). We tested the robustness of our results against this as- 
sumption by repeating the calculations presented in this paper 
with v tur b = or v tur b = v t h- We found no impact on the Lya 
escape fractions. 



3.2. Velocity dispersion of the gas 

The transfer of Lya photons through a parcel of gas depends 
on the velocity distribution of the neutral hydrogen atoms in 
that gas. This dependence appears in the Voigt profile through 
the Doppler parameter b. In the case of purely thermal motion, 
b = v t h = (2&fiT I mm) 112 • In the presence of small-scale turbu- 



3.3. Dust distribution 

In our simulations, gas metallicity is self-consistently calcu- 
lated from the release of metals in the ISM by SN explosion. 
The dust distribution is derived from the metallicity distribution 
by assuming a g alactic value for the dust-to-metal (mass) ratio 
Rdmt/metai = 0.3 (Inoue 2003 ). We further assume that dust is 
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Figure 3. Number of photons emitted per unit velocity bin. The 
black histogram shows the frequency distribution in the sources 
frame, as sampled by all photons used in our G2 experiment. 
This line profile is the one we use as a model for Hn regions 
(Sec. |3.4| ). The dotted red (resp. dashed blue) histograms show 
how these photons are distributed in the observer frame when 
the galaxy is seen face-on (resp. edge-on). 

present only in neutral (i.e. T<10 4 K) media, and write: 



_ Z m H 

^dust — ^dust/metal ~ n Hli 
X H m d 



(3) 



where Z is the metallicity, m H /md = 5 x 10~ 8 is the proton to 
dust particle mass ratio (Dra ine and Lee| fl984), and is the 
numerical density of neutral hydrogen derived in Sec. [5T 



3.4. The Lya sources 

In the present work, the only mechanism which produces Lya 
photons that we take into account is the recombination of photo- 
ionized gas surrounding massive stars ( Kennicutt 1998). In prac- 
tice, given our resolution, we are effectively using young star 
particles as a proxy for Hn regions and shooting photons from 
their locations. Note that we only consider star particles younger 
than 10 Myr as they are the only ones producing significant 
amounts of ionizing photons (e.g. |Charlot and FaH| 1993). 

Each one of these young star particle emits on average 
Afphot/N* photon^] where A/ph t is the total number of photons in 
the experiment, and Af* is the number of young stellar particles 
in the simulation. In our Gl simulation, we fix Af p h t = 6.4 10 5 , 
and since this simulation has Af* ~ 10 3 , this implies an aver- 
age of 640 photons per source. In our G2 simulation, we pick 
Afphot ~ 5.1 10 6 , and since G2 contains Af* ~ 3 10 4 , this yields 
an average of 170 photons per source. 

In the rest- frame of a star particle, the photons are emitted so 
as to sample a profile defined by a flat continuum plus a Gaussian 
line of full- width at half maximum of 20 km s" 1 , and of equiva- 



lent width EW(Lya) = 200 A (as suggested by e.g. Chariot and 
|Fall|1993t [Valls-Gabaud| 19931 |Schaerer|2003|) . The continuum 
extends from -20000 km s" 1 to 20000 km s _i around the Lya 
line in the star particle rest-frame, which is way beyond any res- 
onant transfer effect. As an illustration, we show on Fig. [3] the 

1 We deliberately elect not to fix the number of photons per source, 
but the total number of photons for the whole experiment and sample 
sources randomly. Thus, the number of photons emitted by each indi- 
vidual star particle has Poissonian variance. 



source frame emission profile (black histogram) along with the 
emitted spectrum (emission from all star particles) of the face-on 
(red histogram) and edge-on (blue histogram) G2 galaxy in the 
observer frame. 



4. How the ISM structure impacts Lya transfer 

The different temperature floors of the cooling functions used in 
Gl (10 4 K) and G2 (100 K) result in strikingly different struc- 
turations of the ISM. While the gas in G2 is able to fragment 
into small star-forming clumps, the thermal pressure support in 
Gl yields a rather smooth ISM with homogeneous star forma- 
tion. In the present section, we investigate the influence of these 
structures on the Lya escape fractions and line profiles. 



4.1. Escape fractions 

Let us first look at the escape fraction as a function of emission 
frequency / e sc(v e m), defined as the fraction of photons emitted 
at a frequency v em (in the source rest-frame) which eventually 
escape the galaxy (i.e. which are not absorbed by dust), whatever 
their observed frequency. The black curves on Fig. [4] show / esc 
as a function of v em for Gl (top panel) and G2 (bottom panel). 

Before we discuss these curves in more detail, we can par- 
tially compare our results to those of |Laursen et al.| ( [20 09 ). The 
simulations of these authors are relatively similar to our Gl sim- 
ulation in the sense that their cooling curve also stops at 10 4 K 
( [Sommer-Larsen et al.|2003|) and therefore their ISM is unstruc- 
tured on small scales. As |Laursen et aL] ( |2009| ) only propagate 
Lya photons emitted exactly at the Lya frequency, we have to 
compare their escape fractions to our value at v em = 0, which 
is about 55%. This is broadly consistent with their results (their 
Fig. 9, at M vir = 10 10 Mo). 

We now turn to the comparison of the escape fractions pre- 
dicted for Gl and G2 (black curves in Fig. [4]). The first re- 
markable difference is seen in the continuum escape fractions^] 
which are ~ 95% in Gl versus ~ 22% in G2. This difference 
has a double origin. The first origin is the very different ISM 
structures found in Gl and G2. As gas in G2 is able to cool to 
lower temperatures, it fragments into large star-forming clouds 
within which most young stellar particles are buried. Instead, 
gas in Gl has a stronger pressure support and star formation is 
more diffuse and young stars end up distributed in a rather low 
density environment. Indeed gas density in Gl does not reach 
values above 10 atoms per cm 3 , while gas in G2 reaches 10 4 
atoms per cm 3 . The second origin is that dust, which we model 
to scale with the density of metals (and neutral hydrogen) natu- 
rally follows star formation. Although it mixes on large scales 
with time, the dense star-forming clouds of G2 turn out to be 
more dust-rich than the diffuse (well mixed) ISM of Gl. This 
produces more exti nguished young stellar p opulations in G2, in 
agreement with e.g. Chariot and Fall ( 2000]). 

The second, more subtle point to take form Fig. [4] is that the 
escape fraction of Lya photons is also less in G2 than in Gl, 
even when normalized to that of the continuum. Indeed, the ra- 
tio of line center to continuum escape fractions in Gl and G2 are 
respectively ~ 54% and ~ 20%. We interpret this as a natural ef- 
fect of resonant scattering of Lya photons. In the dense clumps 
of G2 where most photons are emitted, the effect of dust is en- 
hanced by the large Hi column densities which increase the path 



2 Note that the continuum here is only that produced by star particles 
younger than 10 Myr, and not by the full stellar population of each 
galaxy. 
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Figure 4. Escape fraction as a function of emitted frequency 
v em (in the source rest- frame) for simulation Gl (top) and G2 
(bottom). The black curves show the spatially integrated escape 
fractions, i.e. escape fractions averaged over all sources in the 
galaxy. The red (resp. yellow, green) curves show the escape 
fractions for sources in the high (resp. average, low) density re- 
gions. Thumbnails in the lower-left corner of each panel indicate 
the face-on position of sources contributing to each of the three 
curves of the same color. 



of photons and hence their likeliness to hit a dust grain. The Lya 
escape fractions at the line center vary by an order of magnitude 
between our two simulated galaxies, from ~ 50% in Gl to ~ 5% 
in G2. Note that our simulations do not include ionising radia- 
tion from the stellar sources (or from a UV background). This 
missing ingredient could increase the escape fractions, and pos- 
sibly reduce or the difference between Gl and G2. We plan to 
address this issue in a forthcoming paper. 

Fig. [4] also shows color histograms in each panel. These rep- 
resent the escape fraction distributions for different populations 
of stars. The red (resp. yellow, green) curves show the escape 
fractions for photons emitted in high (resp. medium, low) den- 



sity environment^] Thumbnail images in the lower left corner 
of each panel show where the sources contributing to each curve 
are located in the galaxies (see Figs. [T] and [2] for comparison). 
The very strong variation of escape fractions with environment 
in G2, as opposed to the very small variation in Gl, is the key 
to the difference between these two galaxies, and ex plains our 
disagreement with the results of Laursen et al. ( [2009 ). This dif- 
ferent behavior was expected, since the ISM in Gl (along with 
most simulations in the literature) is artificially smoothed on 
small scales by the unrealistically large pressure support asso- 
ciated with 10 4 K gas. Hence, the weak radial variation in the 
escape fraction of Gl simply reflects the gas density gradient 
in that direction. In G2 however, most young stars lie within 



very dense and dusty clouds, as expected from e.g. Chariot and 
Fall| ( |2000| ), and consequently their Lya radiation is heavily ab- 
sorbed. These results demonstrate the necessity to resolve (at 
least some of) the ISM structure of galaxies with hydrodynami- 
cal simulations before post- processing them with Lya radiative 
transfer codes in order to make more realistic predictions for the 
Lya escape fraction. 

It may seem surprising that our results show the opposite 
trend as that expected in the scenario advocated by [Neu feld 
(1991): we find that a clumpy ISM reduces the escape fraction of 
Lya photons relatively to continuum photons instead of enhanc- 
ing it. The reason for this behavior is that our G2 simulation 
predicts a configuration which is quite different from that as- 
sumed by Neuf eld| ( [l991| ). Whereas this author studies the prop- 
agation of photons through a clumpy medium, we find that our 
photons are actually mostly emitted within the clumps, and it is 
their (in)ability to escape from these clumps that dominates the 
results. Clearly, our G2 simulation still suffers from limitations, 
both in terms of the simplified assumptions we make to describe 
the ISM physics (no UV RT but CIE) and in terms of the reso- 
lution, which is insufficient to fully resolve Hn regions. Bearing 
these caveats in mind, we find no support for the 'Neufeld sce- 
nario' . We plan to return to this issue with improved simulations 
in future work. 



4.2. Lya spectra 

We now turn to our predicted spectra, i.e. the distribution of 
observer-frame frequencies of photons which escape Gl and G2. 
We show these spectra in Fig. [5] for three different inclinations: 
face-on (magenta curves), edge-on (blue curves), and seen from 
45 degrees (green curves). Thanks to the symmetry of the prob- 
lem, we checked that the emergent spectra in all azimuthal di- 
rections are the same. The emergent spectra escaping at inclina- 
tions ±cos6 are also identical. To increase the statistics, we then 
sum all photons escaping with \cos6\ in the bin corresponding 
to the chosen inclination, whatever its azimuthal angle is The 
spectra presented on Fig [5] are built by collecting photons within 
cosO bins large enough to a robust signal, but small enough to 
keep all the angular variations (e.g. \cos6\ > 0.95 face-on and 
\cos6\ < 0.05 edge-on). Our spectra are spatially integrated, i.e. 
they are emergent spectra along a given direction but from the 
whole galaxy. 

Looking at the spectra emerging from Gl in the top panel 
of Fig. [5] we see that they are all double-peaked, and roughly 
symmetrical about the line center. The overall intensity (and 
the equivalent width) decreases regularly from the face-on to the 



3 The cuts in local stellar density that we used vary from Gl to G2, 
and are defined arbitrarily to capture the different regimes illustrated in 



Fig- 14] 
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Figure 5. Top : Emergent spectra from Gl along 3 lines of sight 
(6 = in magenta, 6 = n/4 in green, = n/2 in blue). All 
spectra are spatially integrated, i.e. they are emergent spectra 
along a given direction, but from the whole galaxy. Bottom : 
Same as top panel but for G2. 



nent blue peak. Although a very strong large-scale galactic wind 
is present in G2, this outflowing material is mostly ionized and 
tenuous, conditions which render it transparent to Lya photons. 
As can be seen from Fig. [2] (bottom left panel), SN feedback 
also pushes cold neutral gas outside the disc, albeit not very far. 
However, even such a small-scale kinematic feature is enough 
to produce the asymmetry observed in Fig. [5j because the kine- 
matic energy transfered to this neutral gas is significant (cf Fig[7] 
top panel). 

Even if the asymmetry of the G2 spectra qualitatively com- 
pares well to observations (e.g. |Kulas et aL]|2012| Fig 10), the 
shift of the peak and the extension of th e red tail are less im por- 
tant than in most observed spectra (e.g. |Steidel et al.||2010| Fig 
1). We attribute (at least part of) this discrepancy to the mass 
of the galaxy we simulated. This latter is about two orders of 
magnitude less than typical Lyman-break galaxies, and there- 
fore we can hardly expect its kinematics to have the same am- 
plitude. A wrong implementation of the wind mechanism could 
cause the same mismatch. The simplified assumptions used to 
describe the ionisation state of the interstellar gas, i.e. CIE, can 
also lead to an overestimate of the neutral fraction of the inter- 
stellar medium, especially in G2, and affect the spectral shapes. 
However, investingating different feedback recipes and follow- 
ing UV ionising radiation transfer is beyond the scope of this 
paper. 



5. Orientation effects 

We discussed above the impact of the small-scale ISM structure 
on the propagation and escape of Lya photons. In the present 
section, we inspect the global effect of inclination on the ob- 
served Lya properties of our simulated galaxy G2. Indeed, ori- 
entation effects we re predicted by |Charlot and Fallal 993 |);|Chen 
and Neufeld ( 199 4|) and pointed out by recent studiesflLaursen 



and Sommer-Larsen||2007[ |Laursen et al. 
2010[|Barnes et al.|2011||Yajima et al.|20li 



2009; |Zheng et al. 



b). 



edge-on view. This is expected since the typical optical depth 
from the sources to the observer varies in the same way. The 
double-peaked symmetrical line resembles that expected from 
transfer through a static homogeneous medium. This is not sur- 
prising, since the ISM of Gl is rather homogeneous with a ve- 
locity field typical of a globally rotating disk with little internal 
motions (cf Fig|6j bottom panel). The somewhat higher intensity 
in the blue peak (at negative velocities) is the result of transfer 
effects through the gaseous halo, which is in part in-falling to- 
wards the galaxy (cf Fig[6j top panel). Indeed, we checked that 
the spectrum of photons emerging face-on and for which the last 
scattering is in the halo (above the disk plane) is more asymmet- 
ric than that of photons emerging directly from the disk. 

The emergent spectrum from G2 (bottom panel of Fig. [5]) 
is also double-peaked edge-on, but clearly asymmetric face-on. 
The enhancement of the red peak compared to the blue peak is 
a signature of Lya diffusing in an outflowing medium, and in- 
deed, gas in G2 is outflowing at much larger velocities (of order 
a few 100 km/s) than in Gl (where velocities are close to zero 
or inf ailing, see Fig [6] and Fig [7] top panels). We checked that 
if we set the velocity field to km s" 1 in each cell of the simu- 
lation, the spectrum of G2 seen in any direction becomes sym- 
metric about the line center (v = km s" 1 ). We also reversed the 
vertical component of the velocity field, v z , and checked that in 
this case, the emergent spectrum is reversed, with a more promi- 



5.1. Angular escape fraction and angular redistribution 

Let us start by noting that the observed Lya properties of G2 are, 
as expected, invariant around the axis of rotation of the disc, i.e. 
they do not depend on the azimuthal angle. They do however 
strongly depend on the inclination angle 6, which we define here 
as the angle from the rotation axis of the disc (0 = face-on and 
6 = n /2 edge-on). 

The grey curves in Fig. [8] show the angular escape fractions 
of continuum (dashed line) and line (solid line) photon^] These 
escape fractions are computed as a function of the inclination 
at which photons are emitted. The flat solid grey line thus tells 
us that the escape fraction of line photons is isotropic: what- 
ever direction they are emitted along, only about 5% will escape. 
On the contrary, continuum photons tend to escape more easily 
when emitted perpendicular to the disc. Their escape fraction 
varies from 17% when emitted in the plane to 26% when emit- 
ted face-on. The difference between these behaviors is due to 
resonant scattering. Most line photons scatter a huge number 
of times (70% scatter more than 10 6 times, less than 10% scat- 
ter less than 10 times), enough to forget their initial direction. 
On the contrary, continuum photons which escape do so without 



4 We define line photons as having |v em | < 100 km/s in the frame of 
the sources (see black curve of Fig. [3]). For continuum photons, we take 
photons with |v em | > 7000 km/s. 
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Figure 6. Edge-on (top) and face-on (bottom) views of Gl. 
In both panels, the grey scale indicates the density of neu- 
tral Hydrogen (maximum value in a slice of thickness 8 high- 
resolution cells). The arrows indicate the velocity field of the 
gas at the location where photons last scatter before they are ob- 
served (the double sided arrow in the bottom left corner of each 
panel gives the scale), and their color scales with the logio of 
the number of photons seen from each cell (red: lots of photons, 
blue: small number of photons). 



scattering (about 65% escape directly, and all escape with less 
than 8 scatterings). 

The orange curves in Fig. [8] show the distributions of the in- 
clinations with which the photons escape, i.e. the observed incli- 
nation angle distribution, again for continuum (dashed line) and 
line (solid line) photons. The solid orange curve of Fig. [8]shows 
a strong angular redistribution of line photons: although they are 
emitted isotropically and have an isotropic escape fraction, the 
probability (per unit solid angle) of one escaping face-on is about 
15 times higher than that of one escaping edge-on. Similarly, al- 
beit with a much lower amplitude, continuum photons are also 
redistributed in direction, so that they have a probability 3 times 
higher to come out face-on than edge-on. 

To understand better this angular redistribution, we show in 
Fig. [9] the relation between observed (y-axis) and emitted (x- 
axis) inclinations of the photons which escape. The top panel of 
Fig. [9] shows that indeed, most continuum photons are observed 



r a RT in hydro sims 




Figure 7. Same as Fig. |6| for G2. Note the much more pro- 
nounced vertical motions in the top panel. 

with the same inclination as they were emitted with. This panel 
also shows that a small fraction of continuum photons emitted 
edge-on are scattered and observed as coming from any direc- 
tion. The lower panel of Fig. [9] shows that the situation is quite 
different for line photons. Here, photons only have a marginal 
chance to escape in the same direction as they were emitted in, 
provided they were emitted almost face-on, and the overall re- 
distribution is almost independent on the emission inclination. 
This is due, of course, to the very large number of scatterings 
that line photons undergo in general. 

An early prediction of such orientation effect (favoring face- 
on escape) was made by |Charlot and Fall| ( |1993| ), by analogy 
with the stellar 'limb-darkening' effect. However, the enhance- 
ment factor that we get (~ 15) is much higher than their predic- 
tion (~ 2.3). 

This idealised simulation of galaxy formation may be a much 
more symmetrical configuration than galaxy formation in a cos- 
mological context. However, if still relevant for real galaxies, the 
strong orientation effect we find has important observational im- 
plications. In particular, our results suggest that an LAE survey, 
which introduces a hya luminosity selection, will be biased to- 
wards face-on objects, and much more so than e.g. Lyman break 
surveys which rely on a (UV continuum) broad-band selection. 
The LAE may hence represent only an incomplete survey of star 
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Figure 8. Comparison between the distributions of escape direc- 
tions in orange, and emission directions in grey, for line (solid 
lines) and continuum (dashed lines) escaping photons. See the 
text for a detailed explanation of the shapes of these distribu- 
tions. Note that the y axis has two different meanings, f esc for 
the distribution of emission directions, and P for the normalised 
distribution of escape directions. 



forming objects at high redshifts. Also, we expect inclination to 
introduce a significant scatter in correlations between e.g. SFR 
and observed Lya luminosity. 



5.2. Lya equivalent widths 

The inclination effect discussed above also manifests itself in 
the form of a strong correlation between the observed equivalent 
width and the inclination at which the galaxy is observed. In Fig. 



10 



we show the EW as a function of the inclination angle, each 
point corresponding to a random direction of observation. The 
tight correlation is well fit by a polynomial of the form EW(ja) = 
-8 + 80.7 - 393.1 yu 2 + 798.2 ji 3 - 387.4 /A where fi = cos(0). 
This fit is shown with the solid black line. The scatter across 
the relation is due to variations in the azimuthal angle ((f)). Note 
that the EWs we compute here only take into account continuum 
radiation from stars younger than 10 Myr (the same that are used 
as sources for Lya photons), which is good approximation. 

The strong dependency of the EW on \i illustrates the dif- 
ferential effect introduced by resonant scattering and discussed 
with Fig. [8] It shows the complexity of the possible bias in- 
troduced by EW selections in narrow-bancj^LAE samples: not 
only does the Lya luminosity of our simulated galaxy strongly 
vary with inclination, the EW of its emission line does too, and 
as strongly. For example, with a typical selection at EW > 20 A, 
our simulated galaxy would be found as a LAE only if observed 
with an inclination > 0.6, i.e. with ~ 40% chance assum- 
ing a random inclination. In other words, if real galaxies behave 
similarly to our simulation, LAE narrow-band surveys could be 
missing ~ 60% of the faint galaxy population at high redshifts. 
Even if this number is probably an upper limit, given the symetry 
of our idealized galaxy , it calls for observational studies which 




0.0 I 

0.0 0.2 0.4 0.6 
(edge-on) cos(0 em ) 



0.8 1.0 
(face-on) 



5 Although narrow band surveys introduce by construction an EW 
selection, most spectroscopic searches bear such a bias as well. 



Figure 9. Angular redistribution of continuum (top panel) and 
line (bottom panel) photons. The x-axis gives the emission di- 
rection (as the cosine of the inclination angle), and the y-axis 
gives the direction along which the photons escape. The color 
coding scales logarithmically with the number of photon per unit 
area in the plot, as indicated by the color tables. The white lines 
on both panels, show the escape fraction as a function of the 
emission direction. 



could assess the effect of inclination on the Lya properties of 
real galaxies. 



5.3. Line profile 

As presented on Fig [5] the emergent spectrum from our realis- 
tic ISM model G2 shows a transition in its shape with inclina- 
tion: the edge-on profile is a double peak symmetrical around 
the line center, and the face-on spectrum has a much attenuated 
blue peak. The origin of this behavior is mainly illustrated in 
Fig. [7] and is as follows. Lya photons seen escaping edge-on 
have scattered through an extremely thick medium, with little 
radial velocities, and the line they produce thus resembles the 
double peak of a static medium in which photons scatter until 
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Figure 10. Lya equivalent widths as a function of the (cosine 
of the) inclination at which our simulated galaxy is seen. Each 
point marks one of 5000 mock observations. The solid line is a 
simple polynomial fit (see text). 



they are shifted far enough from the line center to escape. Lya 
photons seen face-on, on the contrary, traverse a much lower col- 
umn density with a much more pronounced velocity field, with 
gas outflowing at velocities of several 100 km/s, which favors 
the escape of red-shifted photons. 

This evolution of the line profile suggests that the integrated 
spectral shapes of galaxies may reflect to some extent their in- 
clination, as much as their internal structure and kinematics. 
Clearly, one would need a larger sample of simulated galaxies 
to draw robust conclusions here. However, this trend found in 
our results does hint towards yet another selection effect: Our 
simulation predicts a correlation between the EW and the line 
profile, due to the strong correlation that both these quantities 
have with inclination. This implies that an EW selection will fa- 
vor asymmetric line profiles and tend to exclude double peaks, 
which we find to be associated to low inclinations. 
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6. Lya halo 

On Fig [TTJ wepresent face-on (top panel) and edge-on (bottom 
panel) image^jof the galaxy G2 in the Lya line. Scattering of 
Lya photon through the tenuous intra-halo medium produces a 
diffuse Lya halo which is clearly visible from both angles. The 
face-on view only includes photons with |cos(0)| > 0.95. The 
edge-on view includes all photons with escaping | cos(0)| < 0.2 
so as to have enough statistics, and its surface brightness was 
renormalized accordingly, so that both panels in Fig [TT] are di- 
rectly comparable. The observed Lya luminosity face-on is 
15% of the intrinsic luminosity (cf Fig [8] L(LyoOobs = 0.15 x 
L(Lya){ nt ~ 1.8 x 10 41 erg.s -1 ), whereas it is only a few per- 
cent in the edge-on sightline (~ 1 - 2 x 10 4Q erg.s" 1 ). Su ch faint 
galaxy can be observed from z = (Ostlin et al. 1|2009[ Lyman 
Alpha Reference Sample,), until z ~ 3 (Rauc het al.|2008| [Garel 
et al.|2012| ). As can be seen from the edge-on view, most of the 
extended emission comes from scattering on the envelope of the 



6 Images along a given line of sight described by ki os are obtained 
by selecting photons escaping in a cone of angle a defined by -1 < 
kios'k out < cos(n - or), with a being as small as possible to achieve 
the best accuracy, and high enough to collect a significant number of 
photons. 



Figure 11. Face-on (top) and edge-on (bottom) views of simu- 
lated galaxy G2 in Lya. The color-coding shows surface bright- 
ness in arbitrary units, with the same scale in the two panels. 



galactic wind produced by G2, and the hot wind itself is mostly 
transparent. 

This picture is qualit a tively comparable t o obser vations by 



e.g. 



Hayes et al. (2005); Ostlin et al. (2009) or Steidel et al. 



(201 lb. Quantitatively, we find that about 40% of Lya photons 
come from the central region (within 1.5 kpc), and the rest is 
diffuse. This does not fit well with the results of Ste idel et~aLl 
(2011), who find about 5 times more luminosity in the extended 
halo than in the central region. There are however a number 
of likely explanations to this apparent disagreement. First, our 
simulated galaxy is about two orders of magnitude less massive 
than a typical LBG. Second, and probably more important, our 
simulation is idealized, and the circum-galactic medium does 
not feature cold streams or tidal tails which would enhance the 
effect of scattering. Third, our calculations do not include in- situ 
emission from cold gas in the halo, either from cooling radiation 
or UV background fluorescence, and thi s may play a dominant 
role here (e.g. Rosdahl and Blaizot 2012). 
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7. Summary and conclusion 

In this paper, we have studied the Lya properties of a couple 
of high-resolution simulated dwarf galaxies forming in an ide- 
alized dark matter halo. Our two simulations assume different 
temperature floors of the cooling function (10 4 K for Gl, 100K 
for G2), which result in strikingly different structurations of the 
ISM. While the gas in G2 is able to fragment into small star- 
forming clumps, the thermal pressure support in Gl yields a 
rather smooth ISM with homogeneous star formation. We have 
post-processed these galaxies with McLya in order to follow the 
resonant scattering of Lya photons through their ISM, and to 
predict their resultant observational properties. Our main results 
are as follows. 

1 . As expected, the small scale structure of the ISM plays a de- 
terminant role in shaping a galaxy's Lya properties. In the 
G2 simulation, where gas is allowed to cool down to tem- 
peratures <sc 10 4 K, most young stars are embedded in thick, 
dusty, star forming clouds, and the photons they emit are 
strongly attenuated. As opposed to the "Neufeld scenario", 
the dumpiness of the ISM here enhances the destruction of 
Lya photons relative to continuum photons. This is due to 
th e fact that photons are e mitted within the dense clouds, a 
/fl lCharlot and Fall] ( |2000| ), rather than outside, as assumed 
in [Neufeld| ( |1991| ). 

In the Gl simulation, with an artificially warm ISM, young 
stars are found in much lower density environments and their 
photons escape more easily. This simulation is comparable 
to th e previous studies in the literature which als o include 
dust ( [Laursen et~aL]|2009[ |Yajima et aL]|2011a|b| ), and our 
results are indeed similar to these studies. Such simulations 
do not capture the enhancement of Lya extinction relative to 
continuum in star- forming regions that we find in G2, simply 
because they do not form such dense star forming regions. 
Another important feature is the kinematic structure of the 
ISM. Because G2 develops a genuine multiphase medium, 
with very dense star forming clouds and a tenuous diffuse 
component, supernovae explosions are able to push gas to 
high velocities (see Fig. [7]). Instead, the rather homogeneous 
ISM of Gl is overall denser than the diffuse medium of G2, 
and resists better supernovae explosions. It thus displays a 
rather static velocity field (see Fig. [6]). As shown in Sec. 
|4.2| these different velocity fields have a strong impact on 
the shape of Lya lines. 

2. The analysis of Lya emission from G2^] demonstrates the 
existence of a strong inclination effect. Due to the numer- 
ous scatterings that line photons undergo, their probability 
to escape does not depend on the direction towards which 
they were emitted. Instead, they tend to systematically es- 
cape face-on, following the path of least opacity. Because 
continuum photons do not display such a strong angular re- 
distribution, this effect is directly seen on the Lya equivalent 

width, which we find to vary from 5 A edge-on to ~ 90 

A face-on. 

We also find that this inclination effect is seen in the shape 
of the Lya line emerging from our simulated galaxy. When 
seen edge-on, our galaxy has a double-peak line, associated 
with a low EW. When seen face-on, our galaxy has an en- 
hanced red peak, and a high EW. 

These results suggest the possible existence of strong obser- 
vational biases in LAE surveys which necessarily rely on 
Lya luminosity and EW selections, and could thus prefer- 



entially select face-on objects. As an example, a survey with 
an EW cut at 20 A would select our galaxy only 40% of the 
times, assuming it has a random inclination. 
3. Scattering of galactic Lya photons through the circum- 
galactic medium do produce an extended Lya halo. We find 
that about a third of Lya photons escaping G2 contribute to 
this diffuse component. This is somewhat at odds with the re- 

though the comparison should 



suits of ISteidel et al. (2011 



and Gl, although we do not show it here 



be taken with caution given the fact that our simulated galaxy 
has a much smaller mass than the galaxies analyzed by these 
authors. Also, our simulation is idealized, and the CGM of 
G2 is not a good representation of wha t one finds in cos- 
mological simulations at high redshift ( Rosdahl and Blaizot 
|20l2||Dubois et al.|2012> . 

Although our quantitative results are clearly limited by much 
missing physics, we believe that this work demonstrates the two 
main following points: (1) resolving the ISM is mandatory if 
we want to understand the escape fraction of Lya from galaxies 
(it doesnt matter here if we have the correct solution: we show 
how widely the results vary when we change the ISM physics 
... and this definitely shows that we should go further); and (2) 
for an ideal disc galaxy, we find that the escape fraction is a 
strong function of inclination, and we argue that this effect is 
quite possibly present in real galaxies (eventhough their mor- 
phologies are known to be more complex). Both results call for 
more work, both theoretically and observationally. From the the- 
oretical viewpoint, we plan to make progress in forthcoming pa- 
pers by (i) including the transfer of ionizing photons through a 
resolved ISM, and (ii) embedding our galaxy in the full com- 
plexity of its cosmological context. 
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